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SUMMARY 


The  viscous  transonic  small  disturbance  equation  has  been  recast  into  a  form 
which  displays  the  role  of  the  usual  inviscid  transonic  similarity  parameter  and 
of  a  newly  defined  viscous  transonic  similarity  parameter  in  the  two  dimensional 
flow  past  slender  bodies.  Using  the  approximate  method  of  Hosokawa  the  solu¬ 
tion  for  viscous  transonic  flow  past  a  wavy  wall  has  been  obtained  in  analytical 
form,  and  displays  the  role  of  the  above  similarity  parameters.  The  influence 
of  the  Reynolds  number  and  free  stream  Mach  number  on  the  supersonic  pockets 
which  arise  at  the  wall  surface  has  been  investigated. 
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I.  INTRODUCTION 


Regions  of  supersonic  flow  arise  near  the  surfaces  of  airfoils  and  nozzles 
as  the  free  stream  Mach  number  approaches  unity.  Such  supersonic  pockets 
are  usually,  though  not  always  (Nieuwland  and  Spee,  1968)  terminated  by  a 
weak  shock  wave.  Analysis  of  the  flow  within  and  near  such  pockets,  particu¬ 
larly  in  the  neighborhood  of  the  terminating  shock  is  a  crucial  problem  of  tran¬ 
sonic  aerodynamics.  A  particular  difficulty  is  that  the  pressure  rise  across 
the  terminating  shock  wave  is  not  the  same  as  that  obtained  from  the  Rankine- 
Hugoniot  equations  due  to  the  interaction  between  the  shock  and  the  boundary 
layer  at  the  surface  and  also  due  to  the  fact  that  the  weak  terminating  shock 
may  be  sufficiently  thick  that  the  shock  structure  is  no  longer  one  dimensional 
(Sichel,  1968). 

To  gain  some  insight  into  this  problem  Hosokawa  (1960a,  b),  using  an  ap¬ 
proximate  method,  obtained  solutions  of  the  inviscid  transonic  equation  for 
flow  past  a  wavy  wall.  With  increasing  free  stream  Mach  number  supersonic 
pockets  do  appear  at  the  surface  of  the  wavy  wall,  and  Hosokawa's  analysis 
always  requires  that  such  supersonic  regions  be  terminated  by  a  shock  like 
discontinuity.  Hosokawa's  method  was  applied  by  Sichel  and  Yin  (1969)  to  ob¬ 
tain  solutions  of  the  viscous  transonic  (or  V-T)  equation  for  wavy  wall  flow. 

The  V-T  equation,  which  is  discussed  in  a  survey  by  Sichel  (1968),  corrects 
the  inviscid  transonic  equation  by  an  additional  term  representing  the  influence 
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of  compressive  or  longitudinal  viscosity,  but  shear  viscosity  is  still  neglected. 
Thus  the  V-T  equation  makes  it  possible  to  analyze  weak  shock  waves  with  a 
two  dimensional  internal  structure.  Sichel  and  Yin's  results  could  only  be  ob¬ 
tained  by  numerical  integration  and  hence  were  limited  to  a  single  Reynolds 
number  and  only  several  values  of  free  stream  Mach  number. 

In  the  present  paper  the  V-T  equation  is  first  recast  into  a  form  which 
displays  the  role  of  the  usual  inviscid  transonic  similarity  parameter,  and  a 
newly  defined  viscous  transonic  similarity  parameter  in  transonic  flows  about 
bodies  in  general.  Then,  using  the  method  of  Hosokawa,  the  solution  for  V-T 
flow  past  a  wavy  wall  is  obtained  in  analytical  form  so  that  the  influence  of 
Reynolds  number  and  free  stream  Mach  number  on  the  supersonic  pockets  which 
arise  can  be  explored  in  detail. 
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II.  THE  VISCOUS  TRANSONIC  EQUATION  AND  SIMILARITY 


The  V-T  small  disturbance  equation  has  been  derived  (Sichel,  1968)  from 
the  Navier -Stokes  equations  by  perturbing  with  respect  to  a  uniform  sonic  flow. 
But,  a  more  general  equation  is  obtained  if  perturbations  are  with  respect  to 
some  uniform  free  stream  velocity  U  near  the  sonic  value  but  not  necessarily 
equal  to  it.  Consider  two  dimensional  flow  past  a  body,  whose  surface  is  given 
by 

y-ywW=Sf(x/L)  (1) 


as  shown  in  Fig.  1.  Barred  quantities  are  dimensional.  S  is  the  maximum 

_  _  _  _  _2 

thickness,  and  the  body  shape  function  f(x/L)  ~  0(1).  Then  using  U,  p^,  p^U 
as  reference  velocity,  density  and  pressure;  and  the  body  length,  L,  as  refer¬ 
ence  length,  the  expansions 


u  1  (1) 

u  =—  =  1  +  €  u  + 

U 


=  1  +  e  p 


(1) 


Poo 


+  •  •  • 


+  •  .  • 


(2) 


u  =—  =  Aev 


(1) 


U 


and  the  coordinate  stretching 


x 


L 


(3) 


when  introduced  in  the  Navier  Stokes  equations,  yield  the  following  equations 

,  (1)  ^  (1) 

for  u  and  v  : 
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v"  (.  y  -  l\ 

UlTI  +  Prf'/U 


(1) 


(1  -  *0  (1) 

+ - u 

xx  e  x 


,  1N  ..  2  (1)  (1)  A*  (1)  A 

-  (y  +  1)  M  u  u  +  —  v  =  0 
oo  x  c  y 


(4) 


(5) 


The  stretch  factor  A  takes  into  account  the  difference  in  the  characteristic  x 
and  y  dimensions  which  arises  in  transonic  flow  while  u^,  p^\  and 
are  presumed  to  be  0(1),  and  the  expansion  parameter  e  1.  v"  is  the  com¬ 
pressive  or  longitudinal  viscosity  and  Pr"  is  the  longitudinal  Prandtl  number. 
Since  the  flow  is  irrotational  to  the  present  order  of  approximation  introduc¬ 
tion  of  a  potential  4>  reduces  Eqs.  (4)  and  (5)  to 
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To  the  present  order  of  approximation  the  boundary  conditions  will  be  that 


A£v(\o)=Ae^,0)=^  =  s| 


(7) 


where  S  =S/L,  and  that  u^  and  vu^  must  remain  finite  as  y  —  oo.  At  the  same 


.(1) 


time  the  pressure  coefficient  is  given  by 


4 


ur 


9 


p  -  p 


00 


C  =- - ~  =  -  2c  $ 

p  1-  r,2  x 

2p=oU 


Follov/ing  Ferrari  and  Tricomi's  (1968)  discussion  of  inviscid  transonic  simi¬ 
larity  let 


Ac  =  S 


(8) 


since  both  v^(x,  0)  and  f'(x)  should  be  0(1). 

Using  Eq.  (8)  the  V-T  equations  can  then  be  written  in  the  form 


ULS2  Pr'^ 


(1  -  M  ^  e2 
oo 


4» 


XXX 


(y  +  1)  M  2  c3 
oo 


$  $  +  $  =0 
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(9) 


Now,  in  V-T  flow  the  four  terms  of  Eq.  (9)  should  be  of  :he  same  order,  and  if 
the  expansion  and  stretching  of  Eqs.  (2)  and  (3)  are  appropriate,  the  4>  , 

XXX 

$  ,  $  ,  $  ~  0(1),  and  it  becomes  convenient  to  choose  e  so  that 

xx’  x’  yy 


(y  +  1)  Mw2€3 


= 1  or  e  = 


s2/3 

,  ,.1/3  „  2/3 

(y  +  1) 


(10) 


The  coefficient  of  $  then  becomes 

xx 


(1  -  M  ^e2  (1  -  M  2) 

00  OO 


V/3m4/3(v  +  i)2/3=X“ 


00 
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where  x  is  now  the  usual  inviscid  transonic  similarity  parameter  (Ferrari 
00 


and  Tricomi,  1968).  The  coefficient  of  the  viscous  term  of  Eq.  (9)  becomes 


v"  € 2  f.  r  -  l\ _ v" _  (<  y  -  l)  _  v 

ULS^  Pr’’'~ULS2/3M4/3(,  +  l)2M  PT"> 


(12) 


00 


where  xy  is  now  a  viscous -transonic  similarity  parameter.  The  V-T  equation 


can  thus  be  written  in  the  form 


X  4>  +X$  -  $  4>  +  4>  =0 

v  xxx  oo  xx  x  xx  yy 


(13) 


with  the  boundary  conditions 


$  (x,  0)  sf'(x)  ;  4>  ,  4>  bounded  as  lyl  —  oo,  |x|  -  oo  (14) 

y  y 


The  pressure  coefficient  will  be 


C  = 
P 


2S 


2/3 


,  1/3  2/3  VX,y?Xoo’  V 

(y  +  1)  M 

oo 


(15) 


Equations  (13),  (14),  and  (15),  provide  the  basis  for  viscous -transonic  simi¬ 


larity  rules. 

Since  U  L/v"  is  a  Reynolds  number,  Re,  based  on  compressive  viscosity, 


the  newly  defined  V-T  similarity  parameter,  xy  can  also  be  expressed  in  the 


form 


,2/3  „  4/3 


Re  S  '  M  '  (y  +  1) 
oo 


L  ,  y  - 

2/3  V  Pr"/ 


(16) 
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Thus  xy  varies  inversely  with  the  Reynolds  number,  and  the  coefficient  of  the 
viscous  term  in  Eq.  (13)  will  become  very  small  for  large  Re.  Letting 


n  = 


^273 


u  s 


M 


473 


00 


(y  +  1) 


273 


(••fcO 


(17) 


the  V-T  similarity  parameter  becomes 


xv  -v/h 


(18) 


and  t]  is  a  length  of  the  same  order  as  the  thickness  of  a  weak  shock  wave. 
Finally,  in  terms  of  a  viscous  length  Lv  =  i/’/a^, 

Xv=  3^1737^273  “ 

00 

Since  is  of  the  order  of  the  mean  free  path,  the  ratio  Lv/L  is  clearly  a 
Knudsen  number. 

It  is  readily  shown  that  at  the  sonic  point 


u 


(1) 


d  -  M  J) 

M  4/3  (v  +  l)2/3  '  X=° 
00 


(20) 


that  is  the  critical  value  of  the  expansion  coefficient  u^  equals  the  inviscid 
transonic  similarity  parameter. 
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The  variation  of  x  with  S  and  M  is  shown  in  Fig.  2(a),  while  Fig.  2(b) 

oo  oo 

shows  the  variation  of  xy  with  the  Reynolds  number  U  L /v"  for  different  values 
of  thickness  ratio  S,  and  also  shows  the  variation  of  Lv/L  with  Re.  Clearly, 
the  use  of  the  continuum  theory  is  questionable  when  L*v/L  >  0. 1  or  when 
Re  <  10  in  the  present  case  with  =  1.0.  With  Re  >  100,  Xy  «  1,  except 
for  extremely  small  values  of  the  thickness  ratio  S,  and  then  solution  of  the 
V-T  equation  (13)  becomes  a  singular  perturbation  problem. 

Equation  (15)  implies  the  similarity  rule  that  the  pressure  coefficient  C  , 

2/3  2 

varies  as  S  '  for  fixed  x  ,  X  ,  at  least  to  first  order  in  (1  -  M  )  or  e.  For 

oo  V  oo 

geometrically  similar  bodies  it  is  necessary  to  know  the  influence  of  the  phys¬ 
ical  variables  L,  M^,  S,  and  Lv  upon  the  similarity  parameters  X^and  Xy  to 
assess  the  significance  of  the  V-T  similarity  rule  implied  by  Eqs.  (13),  (14), 
and  (15);  hence  the  effects  of  these  variations  have  been  summarized  in  Table  I 
below.  Case  1  in  Table  I  represents  flow  with  constant  free  stream  density, 
temperature,  and  velocity,  but  with  variable  characteristic  length  L,  which 
might,  for  instance,  be  the  airfoil  chord.  Then  even  though  the  inviscid  tran¬ 
sonic  similarity  parameter  remains  constant,  the  viscous  parameter  Xy  varies 
inversely  with  L.  Flow  past  a  fixed  object  with  constant  p^,  T^  but  variable 
is  represented  by  Case  2,  and  then  xy  is  fixed  but  the  inviscid  parameter 
X^  is  variable.  Cases  3  and  4  show  the  ef’?cts  of  varying  body  thickness  S, 
and  of  varying  p^  and  T^  and  hence  the  viscous  length  upon  the  similarity 
parameters.  Case  5  shows  the  variation  of  C^  when  parameters  are  adjusted  to 

keep  x^  and  fixed. 
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Table  I.  Relation  Between  Physical  Variables  and  Similarity 
Parameters  in  V-T  Flow. 


— 

Case 

No. 

Physical  Variables 

Similarity  Parameters 

Fixed 

Variable 

Fixed 

Variable 

1 

1. 

M  .  S,  L 
ad  ’  v 

L 

*«o 

J^ocL" 

o 

2. 

S,  L,L 

V 

M 

00 

*v 

*00*  (1  -  M»  > 

3. 

L .  i_j  .  M 

’  V1  00 

S 

XcA 

s-rV’ 

4. 

S,  L,  M 

00 

L 

V 

x» 

\ocLv 

5. 

L 

V 

S,  L,M 

00 

l-M^S2/3 

00 

L  oc  S"2 

1 

C  ocS2//3 

P 
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m.  WAVY  WALL  PROBLEM  — FORMULATION 


The  problem  of  flow  past  a  wavy  wall  provides  a  means  of  examining  the 
effects  of  the  parameters  x^  and  Xy  on  V-T  flow  since,  using  the  method  of 
Hosokawa  (I960),  it  is  possible  to  obtain  approximate  solutions  of  the  V-T 
equation  in  analytical  form. 

With  the  ordinate  of  the  wall  given  by 

y  =  S  sin  2ir  ~  (21) 

w  £ 

where  £  is  the  wavelength,  the  boundary  condition  (Eq.  (14))  becomes 

v^  =  $y(x,  0)  =  u>  cos  u>x  (22) 

where  u>=  2rr  L/£.  The  factor  determines  the  choice  of  the  characteristic 
length  L.  Thus  letting  w  =1  implies  that  L  =  £/27t,  and,  as  will  become  evi¬ 
dent  later,  is  the  most  convenient  choice.  In  the  present  formulation  the 

influence  of  S  and  £  enters  only  through  the  similarity  parameters  x  and  x 

oo  v 

in  Eq.  (13),  with  the  choice  of  a  matter  of  computational  convenience,-  in 
contrast  to  the  previous  formulation  (Sichel  and  Yin,  1969)  where  the  param¬ 
eters  of  the  problem  occur  in  both  the  equation  and  the  boundary  conditions. 
Following  Hosokawa  (1960b)  the  potential  is  split  into  two  parts 

*  =  *  +  g  (23) 

such  that  <f>  satisfies  a  linearized  V-T  equation  with  the  acceleration  <f>  in  the 

xx 
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nonlinear  term  replaced  by  a  constant  K  so  that 


X  <f>  +X<£  -  K  <f>  +<f>  =0 

v  xxx  oo  xx  x  yy 


(24) 


<(>  is  made  to  satisfy  the  wavy  wall  boundary  condition  (Eq.  (22))  so  that 


<£y(x,  0)  =  w  cos  cox  =  Re  (co  elwx) 


(25) 


The  function  g  is  a  nonlinear  correction  which  in  view  of  Eqs.  (13),  (24),  and 
(25)  must  satisfy 


Xg  +Xg  +  (K-A)6-(g<b)-gg  +g  =0  (26) 

v  Bxxx  oo  Bxx  yxx  yx  Bx  yx'  Bx  &xx  Byy 

x  33 


gy(x,  0)  =  0 


A  key  assumption  in  the  approximate  analysis,  which  is  discussed  in  detail 
by  Hosokawa  (1960b)  and  by  Sichel  and  Yin  (1969)  is  that 


gyy(x,  0)  =  0  (27) 

so  that  Eq.  (26)  can  be  treated  as  an  ordinary  differential  equation  for  g  at  the 
wall,  y  =0.  Following  Hosokawa  (1960b)  the  constant  K  is  chosen  as  (x,  0) 

XX 

at  the  accelerating  sonic  point,  ^  (x,  0)  =  x  ,  in  the  linear  solution  where 
K  >  0.  This  procedure  yields  the  equations 


<PX  (x*  0; 

<f>  (x*,  0; 

xx  ’  ’ 


X00’XV 


XC0’X 


=  K 


(28) 
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for  K  and  the  sonic  point  x  =  x*.  The  justification  for  the  choice  of  K  is  dis¬ 
cussed  by  both  Hosokawa  (1960b)  and  Sichel  and  Yin  (1969).  Equations  (23)- 
(28)  present  the  basic  elements  of  Hosokawa's  approximate  method. 
Assuming  that  <f>  is  of  the  form 

<f>  =Re  [  F(y)  ela)X] 

the  solution  of  Eq.  (24)  with  boundary  condition  (Eq.  (25))  is  readily  shown  to 
be 

$  =Re £-  (a>/m*^)  exp m1^  y  cos  |)exp  1  (^x  '  y  sin  2  "  f  )  J  (29) 


where 


3  =arctan 


Cl) 


+  wK 

“2 

Cl) 


0  <3  <  7T 


(30) 


Letting  K  =  Xy  =  0  reduces  the  linearized  equation  (21)  to  that  for  inviscid 
linearized  subsonic  or  supersonic  flow,  and  the  solution  (29)  reduces  to  the 
well  known  linearized  subsonic  or  supersonic  wavy  wall  solution. 

Equation  (26)  for  g(x,  0)  can  be  integrated  once  to  yield 

Xvgxx  +  x®Kx'5(,fx  +  gx)2  +  K*=Al  (31) 
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where  Aj  is  a  constant  of  integration.  Introduction  of  the  variable 

C(x)  =gx(x,  0)  +  <£x(x,  0)  -  =  u^(x,  0)  -  (32) 

for  the  deviation  of  u^(x,  0),  the  velocity  at  the  wall,  from  the  critical  value 
and  use  of  the  linearized  solution  for  <f>  then  yields  the  following  ordinary  dif¬ 
ferential  equation  for 


19  A,  1/2 

p  X  u>  ^  oo  X  m 

=  _2^+^+'v 


cos  (ux  +  |  -  |  ) 


(33) 


Equation  (33)  ia  a  Riccati  equation,  which  with  the  transformation 


.0)  3  7T 

«=2X+4+4  ; 


-  1  dT 

-  uxv'T 


changes  to  the  following  second  order  linear  equation  for  T: 


(34) 


T"  +  &  -  2q  cos  20  T  =  0 


(35) 


with 


2A,  -  x 

A  1  00 

a  = 


"2  \2 


q  = 


m 


1/2 


<A  2 


Equation  (35)  is  the  Mathieu  equation  whose  properties  are  well  known 
(McLachlan,  1947).  The  V-T  wavy  wall  flow  has  thus  been  reduced  to  the 
problem  of  solving  the  Mathieu  equation  and  is  considered  below. 
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IV.  RELATION  BETWEEN  X  ,  X  ,  AND  THE  PARAMETERS 

v  oo 

OF  THE  MATHIEU  EQUATION 


The  wavy  wall  solution  depends  upon  the  similarity  parameters  Xy  and 
X^  through  their  influence  on  &  and  q  in  the  Mathieu  equation  (35).  Therefore 
it  is  necessary  to  solve  Eq.  (28),  which  relates  m  and  the  angle  /3  to  Xy  and 

v 

Introducing  the  solution  for  <f>,  Eq.  (28)  becomes 


(36) 


In  solving  Eq.  (36)  it  is  more  convenient  to  deal  with  the  variable 

wx  +  (K/w) 

t\  —  tan  /3  = - - -  (37) 

A 


so  that 

m1/2  =  u>  IxJ^2  (1  +  V2)  (38) 

Upon  squaring  both  sides  Eq.  (36)  can  be  reduced  to  the  following  algebraic 
equation  for  rj 
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Although  Eq.  (39)  is  of  sixth  order  in  rj  so  that  an  analytical  solution  is  not  avail¬ 
able,  a  graphical  interpretation  of  Eq.  (39)  does  provide  an  insight  into  the 
behavior  of  tj.  Let 


(—r) 


then  solutions  of  Eq.  (39)  are  the  intersections  of  the  curves  y^(Tj)  and  ygO?) 
in  the  y-77  plane  as  shown  in  Fig.  3.  As  long  as  w  /\  x  I  >  1.  0,  Eq.  (39) 
will  have  two  solutions  tj^  and  rj2  such  that  <  ^X  /x^  ^  and  from  Eq. 
(37)  it  then  follows  that  K  <  0  when  tj  =  and  K  >  0  when  tj  =  tj^.  Since  K  is 
to  be  evaluated  at  the  accelerating  sonic  point,  tj^  must  be  the  appropriate 
solution  of  Eq.  (39). 

Although  Eq.  (39)  cannot  be  solved  for  tj  as  a  function  of  Xy  and  x^  the 
equation  can  be  solved  for  Xy  and  a  function  of  tj  and  x^  so  that 


and  the  minus  sign  must  be  chosen  in  Eq.  (40)  to  ensure  that  K>  0.  From 

4  16  1/2 

Eq.  (40)  it  follows  that  rj  <  [(<*>  /  1x^1  )  -  1]  if  Xy  is  to  be  real  while 
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j  4/3  2  ^  “ I 

1 7] I  >  |  [(a)  '  /lx  I  )  -  1]  I  if  X  is  to  be  positive.  For  purely  subsonic 
flows  such  that  4>x(x,  0)  <  for  all  x  the  present  formulation  loses  its  mean¬ 
ing. 

The  variation  of  T]  with  xv  and  X^as  determined  from  Eq.  (40)  is  shown 

in  Fig.  4,  for  subsonic  free  stream  flow  with  X  >0,  and  with  w  =  1.  0  for 

00 

convenience  in  calculation.  Figure  4  reflects  the  influence  of  compressive 
viscosity  on  the  linear  solution  for  <f>.  From  the  definitions  of  r\  (Eq.  (37))  and 
q  (Eq.  (35))  it  follows  that 


1 1/2 


q  = 


00 


(i  +  v2) 


1/4 


UX. 


(41) 


and  the  variation  of  q  with  Xy  and  is  shown  in  Fig.  5.  Remarkably,  q  varies 
only  slightly  with  x^  but  is  mainly  a  function  of  the  viscous  parameter  Xy-  For 
a  sonic  free  stream  with  x^  =0,  Eq.  (36)  reduces  to  a  cubic  equation  for  K 
with  solution 


while  the  Eq.  (41)  for  q  becomes 


(Xv  +  K)1/2  1+§XV 

q  = - 2 - = - 2 —  when  xv  ^  1 


(42) 


16 


From  Eq.  (36)  and  (37)  it  can  be  shown  that 


CO 


-1  -1  00 

_  tan  n  +  tan  - 

2  Xi?  -  cox. 


00 


I 


and  the  variation  of  the  accelerating  sonic  point  x*  with  xy  and  x^ 

in  Fig.  6  for  a  subsonic  free  stream  with  x  >0. 

00 


(43) 

is  shown 
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V.  SOLUTION  OF  THE  MATHIEU  EQUATION 


From  Eq.  (41)  for  q  and  the  definition  of  x  and  x  it  can  be  shown  that 

V  00 

q  —  0  as  S  —  0,  that  is  as  the  amplitude  of  the  wavy  wall  decreases  to  zero. 
Solutions  of  the  Mathieu  Eq.  (35)  in  this  limiting  case  q  =  0  will  be  considered 
first  for  guidance  to  the  appropriate  choice  of  the  constants  of  integration. 
Then 


T  =A2  e^  +  A3  e"^ 


where  and  Ag  are  constants  of  integration  and 

\/*»2  -  2Ai ' 


M-  ft- 


«X. 


(44) 


From  Eq.  (34)  it  then  follows  that 


-  »  X,  u  [  (A,/A,)  e"1  -  e’4] 


;=u<\o)-x  - 

(A2/A3)  e*4  - 


(45) 


and  Ag/Ag  is  now  the  only  independent  constant  of  integration.  Letting 
A2/A3  '  “ 


<Z  =  -UXVH 


=  -  -  2A 

»  00  1 


and  choosing  A^  =0  yields 


(46) 
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corresponding  to  undisturbed  free  stream  flow.  This  is  the  expected  result  in 
the  limiting  case  of  vanishing  amplitude.  However,  with  =  0  three  other 
possible  solutions  are 


(a2/a3)  =1.0 

;  ?  -  \o  tanh  u  Xv  * 

(47) 

^A2//A3^  =  ® 

(48) 

(A2/A3)  =-  1.0 

.  th  Xoo  , 

;  C=-X»COtha,x  5 

(49) 

The  solution  (47)  corresponds  to  Taylor’s  (1910)  solution  for  the  structure  of 
a  weak  normal  shock  wave  while  the  solution  (49)  diverges  at  $  =  0.  None  of 
the  solutions  (47) -(49)  are  appropriate  limiting  solutions  for  vanishing  ampli¬ 
tude. 

When  q  /  0,  it  follows  from  the  Floquet  theory  for  second  order  linear 
equations  with  periodic  coefficients  that  solutions  of  the  Mathieu  equation  can 
be  written  in  the  form  (Abramowitz  and  Stegun,  1964) 

T  =A2  e"^  P(£)  +  A3  e ^  P(-£)  (50) 

P(£)  is  a  function  with  period  i r  or  2n  and  n  =  /i(&,q)  is  the  characteristic  ex¬ 
ponent  which  may  be  real,  zero,  imaginary,  or  complex,  and  is  a  function  of 
the  parameters  &  and  q.  With  \i  =  0  one  of  the  solutions  T  will  be  periodic  of 
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period  n  or  2n  and  for  a  given  q  the  values  of  a  resulting  in  /j  =  0  are  called 
characteristic  values  of  the  Mathieu  equation.  With  fi  real  T  diverges  at  either 
£  —  +  oo  or  i  -  -  oo  and  the  solution  is  called  unstable,  while  if  \i  is  imaginary 
T  will  be  periodic  but  of  period  other  than  n  or  2n  and  the  solution  is  termed 
a  stable  solution.  The  (&,  q)  plane  can  thus  be  subdivided  into  regions  of  stable 
and  unstable  solutions  separated  by  fi(&,  q)  =  0.  The  portion  of  the  (&,  q)  plane 
pertinent  here  is  shown  qualitatively  in  Fig.  7  which  also  shows  curves  of 
constant  jx. 

With  q  =  0  it  has  already  been  shown  that  ju  -/T-  X oo/uXy  and  that 
&  <  0  so  that  the  solution  for  T  lies  in  the  unstable  region  along  the  negative 
&  axis  (Fig.  7).  When  q  >  0  it  is  to  be  expected  that  the  solution  will,  at  least 
for  some  range  of  q,  still  lie  in  this  unstable  region,  however  /i  £  V/-^but 
will  depend  on  both  q  and  Then  the  solution  (50)  can  be  expressed  in  the 
form  (McLachlan,  1947) 


rji  A  _/4  n  2rH  ,  *  V  ~  Q 

r  =A2 6  C2re  +  A3e  C2re 

r=-oo  r=-oo 


■2r£i 


with  the  recurrence  relation 


(&-(2r-ifi)  )  C2r -q(C2r+2  +  C2r_2)=0 
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From  Eq.  (52)  it  follows  that  C 2r  and  C _2r  are  complex  conjugates.  Then 
letting 


C 


2r 


Equations  (34)  and  (51)  yield  the  following  solution  for  the  deviation  of  the 
velocity  at  the  wall  from  the  sonic  value: 


"Jv  A3 


00 

+  1 2rp2r  sin  (2r|  +  4,^)  -  pp2r  cos  (2r|  +  <f>2r )  . 

r=l  ^  ‘ 


+  e 


00 

Xf’ 


+  «>o  +  ^  (2rp2r  sin  <2r«  -  4,^)  +  w>2r  cos  (2r«  -  *2j.) 
r=l 


(53) 

)] 


r  00 

i 

0°  ”1 

A2eM 

A, 

p0  +  2  p2r  C0S  (2r?  +  *2r> 

+  e'^ 

J 

p0  +  .5,  p2r  C0S  <2r?  •  V*  I 

3 

L  r=l  J 

i 

r=l  J 

Letting  Ag/A^  -  oo  as  in  the  limiting  q  =  0  case  considered  above  the  solution 
(53)  for  C  becomes 


oo 

2rp2r  3in  <2r^  +  *2^ 

C  =  -  "V  + - ^ - (54) 

p0  +  2,  P2r  008  (2r?  +  V 

r=l 
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From  the  recurrence  relation  (52)  it  follows  that  as  q  -  0,  Cg  and  hence 

p0  -  0  so  that  S  -*  -  o>x  fJi  which  is  identical  to  the  limiting  solution  (46). 
zr  v 

The  solution  (54)  has  period  n  in  i;  or  is  of  period  2ir  in  wx,  that  is  the  same 
period  as  the  wavy  wall. 

If  Ag/Ag  =0  the  solution  reduces  to  Eq.  (48)  in  the  limit  q  -  0.  When 
Ag/Ag  =1.0  the  solution  for 
origin  with  a  superimposed  oscillation,  while  the  solution  will  diverge  at  some 
value  of  4  when  Ag/A^  <  0.  Thus,  the  choice  Ag/A^  -  °o  appears  to  be  the 
only  proper  one  for  the  wavy  wall  problem. 

The  constant  of  integration  determines  the  parameter  ft  (Eq.  (35))  and 
hence  influences  the  solution  through  both  the  recurrence  relation  (52)  and 
through  the  dependence  of  the  characteristic  exponent  p  upon  ft.  With  q  ?  0 
the  choice  of  is  thus  not  as  straightforward  as  in  the  limiting  q  =  0  case. 

A  property  of  Hosokawa's  (1960)  inviscid  analysis  but  not  of  the  present  viscous 
solution  is  that  the  correction  g  vanishes  at  the  accelerating  sonic  point  de- 

X 

fined  by  Eq.  (28).  In  the  viscous  case  g  (x*)  will  depend  upon  the  integration 
constant  A^;  therefore,  in  order  to  compare  the  viscous  and  inviscid  results 
A1  has  been  chosen  to  make  g  (x*)  =0.  Therefore,  at  the  accelerating  sonic 
point,  x  =x*  or  £  =  £*,  where  ?(§*)  =0,the  condition 


behaves  like  a  weak  normal  shock  wave  at  the 


00 


^  2rp2r  sin  (A*  +  *2r) 

p(£,  q)  =  -  (55) 

P0  +  2L  »2r  cos  (2r|*  +  V 
r=l 

must  be  satisfied.  Equation  (55)  is  an  implicit  equation  for  &and  hence  A^, 
which  can  only  be  solved  by  trial  and  error. 

The  series  in  the  solution  (54)  converge  very  rapidly  so  that  only  three 
or  four  terms  of  the  series  need  to  be  retained  in  the  numerical  evaluation  of 
C.  For  low  values  of  q  the  values  of  n  for  different  values  of  &  and  q  have  been 
obtained  from  graphs  in  Abramowitz  and  Stegun  (1964).  For  large  values  of  q 
asymptotic  expressions  given  by  Erdelyi  et  al  (1955)  have  been  used  to  deter¬ 
mine  n  as  described  in  Appendix  A.  In  practice  computations  are  carried 
out  by  first  choosing  a  value  of  q,  which,  from  Fig.  5,  is  equivalent  to  fixing 
the  viscous  similarity  parameter  xy-  C  is  then  computed  for  arbitrary  values 
of  &,  and  the  position  of  the  sonic  point,  £  *  read  from  these  solutions  can  be 
used  to  determine  X  from  Fig.  6.  The  computation  of  C  can  then  be  repeated 
using  values  of  &  corresponding  to  equal  increments  in  x  . 
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VI.  RESULTS  AND  DISCUSSION 


The  variation  of  C  with  x  has  been  determined  for  x  =  0.  61  and  0.  33 

oo  v 

as  shown  in  Figs.  8  and  9.  Hosokawa's  (1960a)  inviscid  solution  for  several 
different  values  of  x  ,  described  in  Appendix  B,  is  shown  in  Fig.  10.  In  the 
case  of  a  sonic  free  stream  with  x^  =  0  the  variation  of  C  is  plotted  in  Fig.  11 
for  Xy  =0.  61,  0.48,  0.  33,  0.  27,  and  for  xy  =  0  corresponding  to  Hosokawa's 
(1960a)  inviscid  solution. 

From  Figs.  8  and  9  it  can  be  seen  that  in  the  viscous  theory  shock  discon¬ 
tinuities  terminating  regions  of  supersonic  flow  are  replaced  by  smooth  compres¬ 
sions  across  which  the  Rankine -Hugoniot  conditions  are  not  necessarily  satisfied. 
Comparison  of  Figs.  8  and  9  shows  the  influence  of  decreasing  xy  or  increas¬ 
ing  Re  upon  the  solutions.  It  can  be  seen  that  the  transition  to  subsonic  flow 
becomes  steeper  with  decreasing  xy-  With  a  sonic  free  stream  x^  =  0,  xy 
does  not  seem  to  affect  the  location  of  the  compression  wave;  in  fact,  the  inflec¬ 
tion  points  of  the  sonic  compressions  occur  at  x  =  1.  25  tt  which  is  also  the 
location  of  the  corresponding  shock  discontinuity  in  the  inviscid  solution  (Fig. 

10).  With  a  subsonic  free  stream,  x^  >  0,  on  the  other  hand  the  effect  of 
viscosity  is  to  shift  the  compressions  upstream  and  to  shorten  the  supersonic 
region  as  compared  to  the  inviscid  solution  of  Hosokawa  (Fig.  10).  Compari¬ 
son  of  the  viscous  and  inviscid  solutions  shows  that  viscosity  has  a  negligible 
effect  on  the  accelerating  portion  of  the  flow  preceding  the  shock  transition. 


24 


Figure  11  shows  that  the  viscous  solution  approaches  the  inviscid  solution 
as  x  -  0  when  the  free  stream  flow  is  sonic  (x  =0);  however,  as  indicated 

V  00 

above  this  no  longer  appears  to  be  true  when  the  free  stream  is  subsonic 

(x  >0).  This  interesting  result,  which  suggests  that  it  may  be  essential  to 
oo 

include  viscous  effects  in  the  analysis  of  certain  transonic  flows,  bears  further 
investigation.  Since  the  theory  developed  here  rests  on  a  number  of  approxi¬ 
mations  and  assumptions,  the  result  that  the  subsonic  solutions  fail  to  approach 
the  inviscid  solutions  as  Xy  -  0  can  only  be  considered  tentative. 

The  sonic  solutions  have  certain  special  features.  For  these  solutions 
the  characteristic  exponent  fi  -  0  so  that  the  sonic  solutions  correspond  to  the 
characteristic  solutions  of  the  Mathieu  equation.  Further,  withXQO=0,  the 
phase  angle  0,  which  arises  in  the  linear  solution  for  <£  has  a  value  of  ir/2  just 
as  in  the  inviscid  case  while  x*  =  tt/4  and  is  independent  of  x  . 

The  choice  of  the  condition  g  (x*)  =0,  used  here,  will  have  a  particularly 
important  influence  on  the  behavior  of  the  subsonic  solutions.  Setting  g  (x*)  =0 
is  ba. J  in  part  on  the  inference  of  Oswatitsch  (1955)  that  the  linearized  solu¬ 
tion  should  be  valid  near  the  accelerating  sonic  point,  and  of  course  this  is  also 
the  point  at  which  the  constant  K  coincides  with  the  actual  acceleration  <f>  of 

XX 

the  linearized  flow;  however,  there  is  no  rigorous  justification  for  this  choice, 
g  (x*)  could  also,  for  example,  be  chosen  to  make  the  inflection  points  of  the 
compressions  in  the  viscous  solutions  coincide  with  the  position  of  the  inviscid 
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shock  discontinuities.  With  this  choice  all  the  viscous  solutions  would  probably 
approach  the  inviscid  solutions  as  xy  0. 

The  flows  considered  here  correspond  to  low  Reynolds  numbers.  For  ex¬ 
ample,  for  a  body  with  thickness  ratio  S  =  0.  01  the  range  0.  27  <  xy  <  0. 61 
corresponds  to  64 >  Re  >28,  and  0.016  <  (Ly/L)  <  0.035  (Fig.  2(b)).  For  this 
value  of  §  will  be  very  close  to  unity  and  the  range  0  <  X^  <  0.  8  corres¬ 
ponds  to  1  >  >  0.970.  Boundary  layers  have,  of  course,  been  neglected 

in  the  present  analysis,  and  the  Re  is  based  upon  the  wavelength  or  character¬ 
istic  body  length.  In  an  actual  flow  the  behavior  of  the  free  stream  shock  wave 
will  also  be  influenced  by  disturbances  induced  in  the  boundary  layer  by  the 
rapid  free  stream  compression.  Then  an  appropriate  Reynolds  number  to  des¬ 
cribe  the  behavior  of  the  flow  near  the  shock  wave  might  be  more  appropriately 
based  upon  a  length  characteristic  of  the  boundary  layer  thickness  rather  than 
the  length  of  the  body.  It  is  thus  possible  even  when  the  Reynolds  number  based 
on  L  is  extremely  large,  that  the  local  Reynolds  number  governing  the  shock 
behavior  may  be  extremely  small,  that  is  of  the  same  order  as  the  Reynolds 
numbers  considered  here. 

As  compared  to  the  previous  paper  on  wavy  wall  flow  (Sichel  and  Yin  1969) 
it  has  now  been  possible  to  obtain  solutions  in  analytical  form.  The  ratio  of 
wall  amplitude  to  wave  length  and  to  a  viscous  length  appears  here  only  through 
the  parameters  x  and  x  whereas  these  ratios  were  introduced  through  the 

V  0C 
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boundary  conditions  in  the  previous  paper.  It  is  difficult  to  compare  the  present 
results  to  those  obtained  by  numerical  integration  in  the  previous  paper  although 
there  seems  to  be  some  difference.  In  particular,  with  a  sonic  free  stream  the 
location  of  the  inviscid  shock  and  the  viscous  compression  do  not  coincide  as  in 
the  present  paper.  As  pointed  out  in  Section  5  above  there  are  several  solutions 
for  £  which  are  periodic  but  only  one  has  the  proper  limiting  behavior  as  the 
wall  amplitude  S  -  0.  While  it  is  easy  to  choose  the  appropriate  solution  here, 
the  difficulty  of  discriminating  between  the  different  solution  branches  when 
numerical  integration  is  used  may  account  for  the  differences  between  the  pre¬ 
vious  and  present  paper. 

Solutions  have  been  obtained  only  for  sonic  and  subsonic  flow,  i.e.  for 
0,  in  the  present  paper.  However  the  present  analysis  can  be  readily 
used  to  evaluate  solutions  for  supersonic  free  stream  flows  with  X^  <  0,  by 
letting  the  constant  of  integration  Ag/A^  =  0  in  Eq.  (53),  and  by  extending  the 
range  of  the  calculations  for  the  relation  between  X  ,  X  ,  and  0.  In  the  limit 

V  QO 

X^  -  0  the  supersonic  solutions  will  approach  the  same  characteristic  Mathieu 
solution  as  the  subsonic  solutions. 

The  present  results  provide  an  indication  of  the  role  of  x  and  x  in  viscous 

V  00 

transonic  flow.  The  boundary  layer  has  been  neglected  and  since  experimental 
results  for  transonic  wavy  wall  flow  are  not  available  there  is  no  way  of  com¬ 
paring  theory  and  experiment. 
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Figure  1.  Definition  of  the  Coordinate  System. 
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Figure  4. 
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Figure  6. 


Figure  7.  Stable  and  Unstable  Regions  in  the  a  -  q  Plane. 
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Figure  9. 


Figure  10.  The  Variation  of  £  with  x;  Hosokawa's  Inviscid  Solution. 
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APPENDIX  A 


For  q  1,  &  >?>  1  the  asymptotic  expressions  for  the  characteristic  value 
given  by  Erdeiyi  et  al  (1955)  can  be  expressed  in  terms  of  complete  elliptic  inte¬ 
grals.  Since  these  expressions  are  not,  to  the  author's  knowledge,  available 
elsewhere  they  are  presented  below. 

For  &  <  -  2  |q  I  the  expression  for  n  can  be  reduced  to 

4  =  (2/ Jr)  (A-l) 

where  E(  )  is  a  complete  elliptic  integral  of  the  second  kind  (Abramowitz  and 
Stegun,  1964).  On  the  other  hand  for  -  2q  <  &  <  2q 

cosh  pr  =  cos  Ij  cosh  I ^  (A-2) 


where 

Ij  =2  ]/ 2q  +  a'[m^  E(l/m)  -  (m  -  1)  m  K(l/m)] 


=  4q/(2q+  a) 

Ig  =2  \J 2q  -  a'fm*^  E(l/m)  -  (m  -  1)  m  ^  K(l/m)] 


m  =  4q/(2q  -  a) 

In  these  expressions  K(  )  is  a  complete  elliptic  integral  of  the  first  kind. 
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APPENDIX  B 


Hosokawa's  (1960b)  inviscid  solution  when  expressed  In  the  variables  used 


here  becomes 


=  ±  \  x  2  + 

*  00 


2K  [*  -*(x*)]  -  2x^x 


(B-l) 


with  the  positive  sign  when  <f>x  >  X^and  the  negative  sign  when  </>x  <  x^.  In 


the  inviscid  case 


-  1/3  ,  1 

<t>  =  -  co  cos  (u)x  -  2  p) 


„  r/  4/3,  2,  ,,1/2 

K=wx00[(w  /X^  )  -  !] 

*(x*)  =  -  (X^/w)  [(a»4/3/xoo2)  -1]1/2 

0  =  arc  tan  [(u>4//3/x  ^  -  lj1^2 
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